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Abstract 

This is a companion to our lectures GAP and loops, to be delivered at the Workshops 
Loops 2007, Prague, Czech Republic. In the lectures we introduce the GAP [B] package 
LOOPS |15| . describe its capabilities, and explain in detail how to use it. In this paper 
we first outline the philosophy behind the package and its main features, and then we 
focus on three particular computational problems: construction of loop isomorphisms, 
classification of small Frattini Moufang loops of order 64, and the search for loops of 
nilpotency class higher than two with an abelian inner mapping group. 

In particular, this is not a user’s manual for LOOPS, which can be downloaded from 
the distribution website of LOOPS. 


1 Maiu features 

On the one hand, since there is no useful representation theory for quasi¬ 
groups and loops, we have decided to represent quasigroups and loops in 
LOOPS by their Cayley tables, thus effectively limiting the scope of the 
package to quasigroups of order at most 500 or so. (A future project is to 
implement other loop representations, notably by connected transversals in 
groups.) 

On the other hand, to take advantage of the powerful methods for groups 
already present in GAP, most calculations in LOOPS are delegated to the 
permutation groups associated with quasigroups, rather than performed on 
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the level of Cayley tables. For instance, to decide if a loop is simple, we 
check whether its multiplication group is a primitive permutation group. 

To avoid repeated calculations, we store most information obtained for a 
given quasigroup as its attribute. In GAP, there is no syntactical difference 
between calling a method or retrieving an attribute. For instance, when Q 
is a quasigroup then CenterCQ) calculates and stores the center Z{Q) of Q 
when called for the first time, while it retrieves the stored attribute Z{Q) 
when called anytime later. 

Moreover, GAP uses simple deduction process—filters—to obtain addi¬ 
tional information about an object without an explicit user’s request. For 
instance, if LOOPS knows that Q is a left Bol loop that is also commutative, 
the built-in hlter (IsMoufangLoop, IsLeftBolLoop and IsCommutative) 
automatically deduces that Q is a Moufang loop and stores this information 
for Q. This is a powerful tool, since many filters built into LOOPS are deep 
theorems. 

1.1 Creating quasigroups and loops 

A {quasigroup) Cayley table is an n x n Latin square with integral entries 
xi < • • • < Xn. A canonical Cayley table is a Gayley table with xi = 1, ..., 
Xn = n. 

When T is a Gayley table, QuasigroupByCayleyTable (T) creates a 
quasigroup whose Gayley table is the canonical Gayley table obtained from 
T by replacing Xj with i. Should T be normalized —the first row and hrst 
column reads xi, ..., Xn —then LoopByCayleyTable(T) returns the corre¬ 
sponding loop. The Gayley table of a quasigroup Q can be retrieved by 
CayleyTable (Q). 

Throughout this paper, we illustrate the methods of LOOPS by ex¬ 
amples, often without any comments for self-explanatory commands. The 
syntax is that of GAP. 

gap> Q := QuasigroupByCayleyTable([[2,1],[1,2]]); Elements(Q); 

<quasigroup of order 2> 

[ ql, q2 ] 

gap> L := LoopByCayleyTableC[[3,5],[5,3]]); Elements(L); L.2; 

<loop of order 2> 

[ 11 . 12 ] 

12 

gap> CayleyTable(Q); 

[ [ 2 , 1 ], [ 1 , 2 ] ] 

gap> Print(L); 

<loop with multiplication table 

[ [ 1 , 2 ], 

[2, 1 ] ] 

> 
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It is also possible to create quasigroups and loops by reading Cayley ta¬ 
bles from files (with very relaxed conditions on the form of the Cayley table), 
by converting groups to quasigroups, by taking subquasigroups, subloops, 
factor loops, direct products, etc. See the manual for details. 

1.2 Conversions 

Even if a quasigroup happens to have a neutral element, it is not considered 
a loop in LOOPS unless it is declared as a loop. Similarly, a group of GAP 
is not considered a loop. We therefore provide conversions between these 
types of algebras: 


gap> G := Group((l,2,3),(1,2)); AsLoop(G); 

Group([ (1,2,3), (1,2) ]) 

<loop of order 6> 

gap> Q := QuasigroupByCayleyTable([[2,1],[1,2]]); AsLoop(Q); 
<quasigroup of order 2> 

<loop of order 2> 


The neutral element of any loop L in LOOPS is always the first element 
of L, i.e., One(L) = L.l. 

Given a quasigroup Q and elements f, g & Q, the principal loop isotope 
{Q,f,g) of Q is obtained from Q via the isotopism {R~^, LJ^, id), cf. |17[ 
p. 60]. Then {Q,f,g) is a loop with neutral element fg. 

The conversion AsLoopCQ) works as follows, starting with a quasigroup 

Q: 


(i) When Q does not have a neutral element, it is first replaced by the 
principal loop isotope {Q,Q.l,Q.l), thus turning Q into a loop with 
neutral element (Q.l)(Q.l). 

(ii) When Q has a neutral element /c, it is replaced by its isomorphic copy 
via the transposition (1, k). 

1.3 Subquasigroups and subloops 

A new quasigroup Q 2 is frequently obtained as a subquasigroup of an ex¬ 
isting quasigroup Qi. Since all information about Q 2 is already contained 
in the Cayley table of Qi, and since it is often desirable to have access 
to the embedding of Q 2 into Qi, we provide a mechanism in LOOPS for 
maintaining the inclusion of Q 2 and Qi. 

When Qi is a quasigroup and S' is a subset of Qi, Subquasigroup(Qi, 
S) returns the subquasigroup Q 2 of Qi generated by S. At the same time. 
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the attribute Parent (( 52 ) is set to Parent (Qi), hence ultimately pointing 
to the largest quasigroup from which Q 2 has been created. The elements of 
Q 2 and the Cayley table of Q 2 are then calculated relative to the parent of 

Q2- 

gap> L := AsLoop(Group((l,2,3),(1,2))): S := Subloop(L,[3]); 

<loop of order 6> 

<loop of order 2> 

gap> Parent( S ) = L; PosInParent( S ); Elements( S ); 

true 

[1,3] 

[ 11. 13 ] 

gap> HasCayleyTable( S ); CayleyTable( S ); 
false 

[ [ 1. 3 ], [ 3, 1 ] ] 

Note that the Cayley table of a subquasigroup is created only upon 
user’s request. 

1.4 Bijections as permutations on {1,..., n} 

When calculating isomorphisms, isotopisms, or other bijections of quasi¬ 
groups of order n, the result is always returned as a permutation (triple of 
permutations) of {1,... ,n}. Equivalently, the quasigroups in question are 
first replaced by isomorphic copies with canonical Cayley tables, and only 
then the bijections are calculated. It is always possible to reconstruct the 
original bijection using the attribute PosInParent. 

1.5 A few words about the implementation 

One of the biggest strengths of the computer algebra system GAP is that 
most algebraic structures can be defined within it. In this subsection we 
briefly explain how the variety of quasigroups is implemented in LOOPS. 
In order to understand the implementation, we will need the following GAP 
terminology: 

- A filter, such as Isinteger and IsPermGroup, is a special unary func¬ 
tion on the set of GAP objects which returns either true or false. 
Roughly speaking, a hlter is an a priori attribute of an object. 

- A category is a class of objects dehned by a collection of Liters. An 
object can lie in several categories. For example, a row vector lies in 
the categories IsList and IsVector. 
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- All GAP objects are partitioned into families. The family of an object 
determines its relation to other objects. For instance, all permutations 
form a family, and groups presented by generators and relations form 
another family. However, a family is not a collection of objects, but 
abstract information about objects. 

- Beside its name, a family can have further labels. 

- Every GAP object has a type. The type of an object determines if a 
given operation can be performed with that object, and if so, how it 
is to be performed. The type of an object is derived from its family 
and its filters. 

- A given data structure can be made into an objeet by specifying its 
type, that is, its family and its filters. 

The following function constructs a quasigroup Q with Cayley table ct. 
First we define a family corresponding to the elements of Q and tell GAP 
that it will consist of quasigroup elements. Then we objectify the individual 
elements in this family, and label the family by the set of its elements, by the 
size of Q, and by the Cayley table. Then we objectify Q whose family will 
be the collection of its elements. Finally, we set some important attributes 
of Q. 


functionC ct ) 

local F, Q, elms, n; 

# constructing the family of the elements of this quasigroup 

F := NewFamilyC "QuasigroupByCayleyTableFam", IsQuasigroupElement ); 

# installing data ("labels") for the family 
n := Length ( ct ); 

F! .size := n; 

elms := ImmutableC List( [l..n], i -> ObjectifyC 

NewTypeC F, IsQuasigroupElement and IsQuasigroupElmRep), [ i ] ) ) ); 
F!.set := elms; 

F!.cayleyTable := ct; 

# creating the quasigroup by turning it into a GAP object 

# the family of Q is the collection of its elements 
Q := ObjectifyC NewType( FamilyObj( elms ), 

IsQuasigroup aind IsAttributeStoringRep ), rec() ); 

# setting some attributes for the quasigroup 
SetSizeC Q, n ); 

SetAsSSortedList( Q, elms ); 

SetCayleyTableC Q, ct ); 
return Q; 

end; 


Operations in GAP are overloaded, i.e., the same operation can be ap¬ 
plied to different types of objects. In order to deal with this situation, GAP 
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uses a method selection: When an operation is called, GAP first checks the 
types of the arguments, and then selects the appropriate method. 

Here is how the multiplication of two quasigroup elements is imple¬ 
mented: 


InstallMethodC \*, "for two quasigroup elements", 
IsIdenticalObj, 

[ IsQuasigroupElement, IsQuasigroupElement ], 
functionC x, y ) 
local F; 

F := FamilyObj( x ); 

return F!.set[ F!.cayleyTable[ x![ 1 ] ][y![ 1 ] ] ]; 
end ); 


Note that the underlying quasigroup is easily accessed since the element 
X knows into which quasigroup it belongs. 

2 What is in the package 

Here is a very brief overview of the methods implemented in LOOPS, ver¬ 
sion 1.4.0. See the manual for (much) more details. Argument Q stands 
for a quasigroup, and L for a loop. Thus the methods with argument Q 
apply to both quasigroups and loops, while those with argument L apply 
only to loops. Any additional restrictions on the arguments are listed in 
parentheses. The symbol > is a shortcut for returns. 

2.1 Basic methods and attributes 

Cayley tables and elements: 

Elements (Q) > list of elements of Q, 

CayleyTable (Q) > Cayley table of Q, 

One(L) 0 the neutral element of L, 

MultiplicativeNeutralElement (Q) > the neutral element of Q, or fail 
SizeCQ) > the size of Q, 

Exponent (L) > the exponent of L {L power-associative). 

Arithmetic operations: 

LeftDivisionCx,y) > x\y, 

RightDivisionCx,y) > xfy, 

LeftDivisionCayleyTableCQ) > Cayley table of left division in Q, 
RightDivisionCayleyTableCQ) > Cayley table of right division in Q. 
Powers and inverses: 

Leftlnverse(x) l> x^, where x^x = 1, 
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Rightinverse (x) l> x^, where xx^ = 1, 

Inverse (x) > the two-sided inverse of x, if it exists. 

Associators and commutators: 

Associator (x, y, z) l> the unique element u with {xy)z = {x{yz))u, 
Commutator (x, y) l> the unique element v with xy = {yx)v. 

Generators: 

GeneratorsOf Quasigroup(Q) l> a generating subset of Q, 
GeneratorsOfLoopCD > a generating subset of L, 

GeneratorsSmallest (Q) > a generating subset of Q of size < log 2 \Q\- 
Subquasigroups: 

IsSubquasigroupCQ, S) > true if 5 is a subquasigroup of Q, 
IsSubloopCL, S) > true if 5 is a subloop of L, 

AllSubloops (L) 0 list of all subloops of L, 

RightCosets (L, 5) > right cosets modulo S {S < L), 

RightTransversal (L, 5) l> a right transversal modulo S {S < L). 
Translations and sections: 

LeftTranslationCQ, x) l> the left translation Lx by x in Q (x G Q), 
RightTranslationCQ, x) > the right translation Rx by x in Q (x G Q), 
LeftSectionCQ) > the set of all left translations in Q, 
RightSectionCQ) > the set of all right translations in Q. 

Multiplication groups: 

LeftMultiplicationGroupCQ) > the left multiplication group of Q, 
RightMultiplicationGroupCQ) > the right multiplication group of Q, 
MultiplicationGroupCQ) > the multiplication group of Q, 
RelativeLeftMultiplicationGroupCL, 5) l> the group generated by all 
left translations of L restricted to S' {S < L), 
RelativeRightMultiplicationGroupCL, S) > the group generated by all 
right translations of L restricted to S {S < L), 
RelativeMultiplicationGroupCL, S) > the group generated by all 
translations of L restricted to S {S < L). 

Inner mapping groups: 

InnerMappingGroup(L) > the inner mapping group of L, 
LeftlnnerMappingGroup(L) > the group generated by L~xLyLx, 
RightInnerMappingGroupCD > the group generated by R~y RyRx- 
Nuclei: 

LeftNucleusCQ) > the left nucleus of Q, 

RightNucleus (Q) > the right nucleus of Q, 

MiddleNucleus (Q) > the middle nucleus of Q, 

NucCQ), NucleusQf Quasigroup(Q) > the nucleus of Q. 
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Commutant, center and associator subloop: 

Commutant (Q) > {x G Q] xy = yx for every y G Q}, 

Center (Q) l> the center of Q, 

AssociatorSubloop(L) > the smallest S < L such that L/S is a group. 
Normal subloops: 

IsNormaKL, S) > true if 5 is a normal subloop of L, 

NormalClosure(L, S') l> the smallest normal subloop of L containing S, 
IsSiinple(L) > true if L is a simple loop. 

Factor loops: 

FactorLoopCL, A^) > L/N (N normal subloop of L), 
NaturalHomomorphismByNormalSubloopCL, A^) > the projection 
L —7> L/N (N normal subloop of L). 

Central nilpotency and central series: 

NilpotencyClassQfLoop(L) > the (central) nilpotency class of L, 
IsNilpotent (L) > true if L is nilpotent, 

IsStronglyNilpotent (L) > true if the mult, group of L is nilpotent, 
UpperCentralSeries (D l> the upper central series of L, 
LowerCentralSeries (L) l> the lower central series of L, 

Solvability: 

IsSolvable(D l> true if L is solvable, 

DerivedSubloopCD l> the derived subloop of L, 

DerivedLength(L) > the derived length of L, 

FrattiniSubloop(L) > the Frattini subloop of L (L strongly nilpotent). 
Isomorphisms and automorphisms: 

IsomorphismLoops (L, M) l> an isomorphism of loops L —>■ M, or fail, 
LoopsUpToIsoinorphism(/s) l> filtered list Is of loops up to isomorphism, 
AutoinorphismGroup(L) > the automorphism group of L, 
IsomorphicCopyByPerm((5,p) > an isomorphic copy of Q via the 
permutation p, 

IsomorphicCopyByNormalSubloopCL, S') l> an isomorphic copy of L in 
which S<L occupies the first |S| elements of L and where the remaining 
elements correspond to the cosets of S in L. 

Isotopisms: 

IsotopismLoops (L, M) l> an isotopism L — M, or fail, 
LoopsUpToIsotopism(/s) > filtered list Is of loops up to isotopism. 

2.2 Testing properties of quasigroups and loops 

Associativity, commutativity and generalizations: 

IsAssociativeCQ) > true if Q is associative. 
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IsCommutativeCQ) > true if Q is commutative, 

IsPowerAssociative (D l> true if L is power associative, 
IsDiassociative(L) > true if L is diassociative. 

Inverse properties: 

HasLeftInverseProperty (L) l> true if x^{xy) = y, 
HasRightlnverseProperty(L) > true of {yx)xP = y, 
HasInverseProperty (L) i> true if x^{xy) = y = {yx)x^, 
HasTwosidedInverses (L) > true if x^ = x^, 

HasWeakInverseProperty (L) l> true if [xy)^x = y^, 
HasAutomorphicInverseProperty(L) > true if {xy)^ = x^y^, 
HasAntiautomorphicInversePropertyCD > true if (xy)^ = y^x^. 

Some properties of quasigroups: 

IsSemisymmetricCQ) > true if {xy)x = y, 

IsTotallySymmetric (Q) > true if Q is semisymmetric and commutative, 
Isidempotent (Q) l> true it = x, 

IsSteinerQuasigroupCQ) l> true if Q is totally symm. and commutative, 
IsUnipotent (Q) > true if x^ = 

IsLeftDistributive(Q) l> true if x{yz) = {xy){xz), 
IsRightDistributiveCQ) > true if {xy)z = {xz){yz), 

IsDistributive (Q) > true if Q is left and right distributive, 
IsEntropic(Q), IsMediaKQ) > true if {xy){uv) = {xu){yv). 

Loops of Bol-Moufang type: 

IsExtraLoop(L) > true if x{y{zx)) = {{xy)z)x, 

IsCLoop(L) 0 true if x{y{yz)) = {{xy)y)z, 

IsMoufangLoop(L) > true if {xy){zx) = {x{yz))x, 

IsRCLoopCD > true if x{{yz)z) = {xy){zz), 

IsLCLoopCD > true if {xx){yz) = {x{xy))z, 

IsRightBolLoopCD l> true if x{{yz)y) = {{xy)z)y, 

IsLeftBolLoop(L) > true if x{y{xz)) = {x{yx))z, 

IsFlexible (Q) > true if x{yx) = (xy)x, 

IsRightNuclearSquareLoop(L) l> true if x{y{zz)) = {xy){zz), 
IsMiddleNuclearSquareLoop(L) > true if x{{yy)z) = {x{yy))z, 
IsLeftNuclearSquareLoop(L) > true if {xx){yz) = {{xx)y)z, 
IsRightAlternative (Q) > true if x{yy) = {xy)y, 
IsLeftAlternativeCQ) > true if {xx)y = x{xy), 

IsAlternative(Q) > true if it is both left and right alternative. 

Power alternative loops: 

IsLeftPowerAlternative (L) l> true if x"'(x”^y) = x"'~^'^y, 
IsRightPowerAlternative(L) > true if {xy"')y"^ = 
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IsPowerAlternative (L) l> true if L is left and right power alternative. 
Conjugacy closed loops: 

IsLCCLoop(L) > true if left translations are closed under conjugation, 
IsRCCLoop(L) l> true if right translations are closed under conjugations, 
IsCCLoopCD l> true if L is left and right conjugacy closed. 

Additional varieties of loops: 

IsLeftBruckLoop(L), IsLeftKLoop(L) > true if L is left Bol and has 
the automorphic inverse property, 

IsRightBruckLoopCD, IsRightKLoop(L) > true if L is right Bol and 
has the automorphic inverse property. 

Here is a nice, albeit trivial illustration of the filters built into the LOOPS 
package: 


gap> L := LoopByCayleyTable([[1,2],[2,1]]); 

<loop of order 2> 

gap> IsLeftBolLoop(L); L; 

true 

<left Bol loop of order 2> 
gap> IsRightBolLoop(L); L; 
true 

<Moufang loop of order 2> 
gap> IsAssociative(L); L; 
true 

<associative loop of order 2> 


2.3 Libraries 

Several libraries of small loops up to isomorphism are included in LOOPS. 
As of version 1.4.0, the libraries contain: 

- all nonassociative left Bol loops of order < 16, 

- all nonassociative Moufang loops of order < 64 and = 81, 

- all nonassociative Steiner loops of order < 16, 

- all (three) nonassociative conjugacy closed loops of order for every 
odd prime p, 

- all (one) nonassociative conjugacy closed loops of order 2p, for every 
odd prime p, 

- the smallest nonassociative simple Moufang loop (of order 120), 

- all nonassociative loops of order < 6. 
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There is also a library of all nonassociative loops of order < 6 up to iso- 
topism. 

The mth loop of order n in a given library can be retrieved via 
LeftBolLoop (n, m), Mouf angLoop(n, m ), 


and so on. 

We took great care to store the information in the libraries efficiently. 
For instance, the library of Moufang loops can be packed into less than 18 
kilobytes, hence averaging about 4 bytes per loop. 

Remark 2.1. All nonassociative Moufang loops of order less than 64 can 
be found in m- Our numbering for these loops agrees with jQ/. 

The Jf262 nonassociative Moufang loops of order 64 were first constructed 
in m, but it was proved (computationally) only in |16]/ that the list is 
complete. 

The 2038 nonassociative left Bol loops of order 16 were enumerated for 
the first time by Moorhouse m- The first author obtained the same result 
by a different method, on which he will report in a separate paper m 

The fact that for every odd prime p there are precisely three nonasso¬ 
ciative conjugacy closed loops of order p^ was established by Kunen m- 
Drdpal and Csorgo derived simple formulas for multiplication in these three 
loops When p is an odd prime, Wilson m constructed a nonassociative 
conjugacy closed loop of order 2p, and Kunen m showed there are no other 
such loops. 

Our counts of small loops agree with the known results, e.g. m 

The library of small Steiner loops is based on m- 


3 Constructing isomorphisms 

There does not appear to be much research on the problem of hnding an 
isomorphism between loops. In this section we explain the approach used 
in LOOPS. It works surprisingly well for many varieties of loops, including 
Moufang loops. 

Let Q be a loop, and let V he a set of properties (of elements) invariant 
under isomorphisms. The nature of V depends on Q. For instance, when Q 
is power-associative, one of the invariant properties for an element x might 
be the order |x|. 

Given V and a collection C of loops, dehne an equivalence on the (dis¬ 
joint) union of C by x ~ y if and only if (p{x) = (p{y) for every € P. 
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Then, if / : Q —>■ L is an isomorphism and C = {Q, L}, we must have 
X ~ f{x) for every x £ Q. In other words, V partitions the elements into 
blocks invariant under isomorphism. 

In order to find an isomorphism, we need a set of invariants V that is 
easy to calculate but results in a fine partition. 

We have used the following invariants V for power-associative loops: 

iplix) = |x|, 

^ 2 {x) = \{y; = x}|, 

= \{y; y^ = x}|, 

I>i,k{x) = \{y, xy = yx, \y\ = k}\, for A: > 1. 

The algorithm searching for an isomorphism f : Q ^ L hrst orders the 
equivalence classes of ~ by increasing size on both Q and L. If the equiva¬ 
lence class sizes of Q and L do not match, it is clear that no isomorphism 
f : Q ^ L exists, and fail is returned. Otherwise, a backtrack search 
attempts to hnd an isomorphism respecting the partitions of ~. 

It would be an interesting project to analyze the speed of the algorithm 
depending on the choice of V. We do not claim that the above V is optimized 
in any sense. Note, for instance, that the invariants ip 2 , are useless for 
many power associative loops of odd order, and are useless for all 
commutative loops. 

4 Classification of small Prattini Moufang loops of 
order 64 

Let L be a loop and let the Frattini subloop d’(L) be the normal subloop 
generated by all squares, commutators and associators of L. In other words, 
<h(L) is the smallest normal subloop such that L/d>(L) is an elementary 
abelian p-group. Following Hsu |9], we say that L is a small Frattini p-loop 
if |d>(L)| < p. 

In this section, L will denote a small Frattini Moufang 2-loop of order 
2 >t--i-i_ Moreover, in order to avoid trivialities, we assume that |‘h(L)| = 2. 
Clearly, d>(L) < Z{L), L is nilpotent of class 2, and it has a unique nontrivial 
square, commutator and associator element. 

Remark 4.1. Small Frattini Moufang 2-loops are also called code loops due 
to their connection to doubly even linear binary codes. Some of these loops 
play an important role in the description of large sporadic simple groups. 
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We consider V = L/d>(L) as a vector space over F 2 , and we identify d’(L) 
and F 2 . In particular, we sometimes write the group operations additively. 

Let us take arbitrary elements u = x mod d>(L), v = y mod ‘h(L), 
w = z mod <h(L) of V. Then, the following maps are well defined: 

a : V ^ F2, cj(m) = 

k:VxV^¥2, k{u,v) = [x,y], 

a \V xV xV ^ F2, Q;(ti, u, w) = [x, y, z]. 

Moreover, a is an alternating trilinear form, k is alternating, and we have 

a{u + u) = a{u) + a{v) + k{u, v), 

k{u + v,w) = K,{u, w) + k{v, w) + a{u, v, w). 

Hence, by definition, H is a symplectic cubic space. 

There are different ways in which a small Frattini Moufang 2-loop is 
obtained from a symplectic cubic space (cf. Griess [8], Chein and Goodaire 
m, Hsu El). All of the above constructions induct on the dimension of V. 
In contrast, a new approach, |13j . takes advantage of groups with triality 
and constructs the loop globally. 

For this, let Uj, Kij and aijk be the structure constants of a,K,a with 
respect to a fixed basis ov V. We define the group G with gerenators 
9i: fij hi, i G {1 ,..., n}, u and v by the following relations: 


9‘f = u^\ 

f! 

= , hj = u^ = v 

" = 1, 

\gi ^ gj\ ~ 




[giJj] = 

{uv\ 

71 




k=l 


[gi,hj] = 


, [fi,hj] = 


\hi, hj\ — 

bo 

II 

II 

= [9hv\ 


Then, G is a group and the maps 

T ■■ gi fi, hi hi, u ^ V 

P : Qic^ fi, fi^ {gifi)~^, hic^ hi, V, V uv 

extend to triality automorphisms of G. The following function returns the 
Moufang loop associated to the group G with triality automorphisms r, p: 
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TrialityGroupToLoop := function( G, tau, rho ) 
local ccl, ct; 

ccl := ElementsC CoiijugacyClass( G, tau ) ); 
ct := List( ccl, si -> 

List( ccl, s2 -> 

PositionC ccl, sl''rho * s2''(rho‘"2) * sl''rho ) 

) 

); 

return LoopByCayleyTable( NormalizedQuasigroupTable( ct ) ); 
end; 


To complete the classification of small Frattini Moufang 2-loops of order 
64, it now suffices to classify the symplectic cubic spaces of order 32. For a 
fixed basis, such a space is given by 



structure constants, which give rise to a 25-dimensional vector space W over 
F2. 

Any linear map A of V defines a new symplectic cubic space with maps 

cr^{u) = a{Au), K^{u,v) = k{Au,Av), a^{u) = a{Au, Av, Aw), 

and hence A induces a linear map on W. This defines an action of GL{5, 2) 
on W. 

It is easy to show the 1-1 correspondence of loop isomorphisms and 
linear isomorphisms of symplectic cubic spaces. This implies that the orbits 
of GL(5, 2) on W will correspond precisely to the isomorphism classes of 
small Frattini Moufang 2-loops of order 64. 

Since |GL(5,2)| and 2^® are still too large for GAP to compute the 
needed orbits, one has to have a closer look at invariant subspaces of W. 
Once this is done, the classification is complete, with the result that there 
are precisely 80 nonisomorphic small Frattini Moufang loops of order 64. 


5 An interesting Csorgo loop 

One of the longer-standing problems in loop theory was the question if 
there is a loop with nilpotency class higher than two whose inner mapping 
group is abelian. In [3], Csorgo constructed such a loop (of order 128 and 
nilpotency class 3). The following GAP code returns this loop L. The code 
follows [3], where some insight is given. 
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# constructing a group of order 8192 by presenting relations 
f := FreeGroup(13); 

G := f/[ f.1-2, f.2-2, f.3-2, f.4-2, f.5-2. f.6-2, f.7-2. f.8-2, f.9-2. f.10-2, 
f.11-2. f.12-2. f.13-2. (f.l*f.2)-2, (f.l*f.3)-2. (f.l*f.4)-2, (f.l*f.5)-2, 

(f .l*f .6)-'2, (f.l*f.7)-2, (f.l*f.8)-2. (f .l*f .9)-'2, (f . l*f . 10) “2. (f.l*f.ll)-2. 

(f.l*f.12)-2, (f.l*f.13)-2. (f.2*f.3)-2. (f.2*f.4)-2, (f.3*f.4)-2. (f.2*f.5)-2. 
(f.2*f.6)-2, (f.2*f.7)-2, (f.3*f.5)-2. (f.3*f.6)-2, (f.3*f.7)-2, (f.4*f.5)-2. 
(f.4*f.6)-2, (f.4*f.7)-2, (f.2*f.9)-2. (f.2*f.10)“2. (f.3*f.8)-2. (f.3*f.10)“2, 
(f.4*f.8)-2, (f.4*f.9)-2, f.l*f.2*f.8*f.2*f.8, f.l*f.3*f.9*f.3*f.9, 
f.l*f.4*f.10*f.4*f.10, (f.2*f.ll)-2, (f.2*f.12)-2, (f.2*f.13)*2, (f.3*f.11)-2, 

(f.3*f.12)-2, (f.3*f.13)-2. (f.4*f.ll)-2, (f.4*f.12)*2, (f.4*f.13)*2, (f.5*f.6)-2. 
(f.5*f.7)-2, (f.6*f.7)-2, (f.5*f.9)-2. (f.5*f.10)“2. (f.6*f.8)-2. (f.6*f.10)“2. 

(f .7*f .8)-'2, (f.7*f.9)-2, f .l*f .5*f .8*f .5*f .8. f . l*f . 6*f . 9*f . 6*f . 9, 
f.l*f.7*f.10*f.7*f.10, (f.5*f.12)-2, (f.5*f.13)-2, (f.6*f.11)*2, (f.6*f.13)*2, 

(f.7*f.ll)-2, (f.7*f.12)-2, f.l*f.ll*f.5*f.ll*f.5, f.l*f.12*f.6*f.12*f.6, 
f . l*f . 13*f . 7*f . 13*f . 7, f . 2*f . 5*f . 9*f . 10*f . 9*f . 10, f . 3*f . 6*f . 8*f . 10*f . 8*f . 10, 
f.4*f.7*f.8*f.9*f.8*f.9, (f.8*f.ll)-2, (f.9*f.12)“2, (f.10*f.13)“2, 
f . 8*f . 12*f . 8*f . 4*f . 12*f . 7, f . 8*f . 13*f . 8*f . 3*f . 13*f . 6, f . 10*f . 1 l*f . 10*f . 3*f . 1 l*f . 6, 
f . 9*f . ll*f . 9*f . ll*f . 7, f . 9*f . 13*f . 9*f . 13*f . 5, f . 10*f . 12*f . 10*f . 12*f . 5, 

(f.ll*f.12)-2, (f.ll*f.13)-2, (f.12*f.13)-2 ]; 

# auxiliary data 

g := GeneratorsOfGroup(G); 

N := Subgroup( G, [ g[5] , g[6] , g[7] , g[l] ] ); 

W := Subgroup ( G, [ g[5]*g[2], g[6]*g[3], g[7]*g[4], g[l] ] ); 

A_0 := [ One(G), g[8] , g[9] , g[10] , g[8]*g[9], g[8]*g[10], g[9] *g [10] *g [2] , 
g [8] *g [9] *g [10] *g [2] ] ; 

B_0 := [One(G), g[8]*g[ll], g[9]*g[12], g[10]*g[13], g[8] *g [11] *g [9] *g[12] , 
g[8]*g[ll]*g[10]*g[13]*g[3] , g[9]*g[12]*g[10]*g[13] , 
g[8]*g[ll]*g[9]*g[12]*g[10]*g[13]*g[3] ] ; 

A := Union( List( ElementsC N ), x -> A_0*x ) ); 

B := Union( List( ElementsC W ), x -> B_0*x ) ); 

H := SubgroupC G, [ g[2] , g[3] , g[4] , g[ll], g[12] , g[13] ] ); 

# constructing the loop 
ListPosition := functionC S, x ) 

local i; i := 1; while not x in S[i] do i := i + 1; od; return i; 

end; 

m := MappingByFunctionC Domain(ElementsC G)), DomainC[1..8192]), 

X -> PositionC Elements(G), x ) ); 

CA := ListC A, X -> x*Elements( H ) ); 
mCA := ListC CA, c -> Set( c, x -> x~m ) ); 

T := ListC[1..128],i->[l..128]); 

for ii in [1..128] do for jj in [1..128] do 

T[ii] [jj] := ListPositionC mCA, (A[ii] *B[j j] )“m ); 
od; od; 

p := SortingPermC T[l] ); 

T := ListC T, r -> PermutedC r, p ) ); 

L := LoopByCayleyTableC T ); 


In addition, the following properties hold for L: (a) the nucleus of L is 
elementary abelian of order 16, (b) the left and middle nuclei have order 
32, (c) the right nucleus has order 16, (d) the two-element center coincides 
with the associator subloop. 

An interesting, more symmetric loop K is obtained from L by this greedy 
algorithm: 
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Given a groupoid Q, let = |{(a, b,c) € Q x Q x Q; a{hc) ^ (a6)c}|. 
Hence ^{Q) is a crude measure of (non)associativity of Q. 

Let T be a multiplication table of L split into blocks of size 16 x 16 
according to the cosets of the nucleus of L. Let h be the nontrivial central 
element of L. 

(*) For 1 < i 7 ^ J < 16, let Tij be obtained from T by multiplying the 
(i,j)th block and the {j,i)th block of T by h. Let {s,t) be such that ^-{Tst) 
is minimal among all fj,{Tij). If ^(Tst) > stop, and return T. Else 

replace T by Tgt, and repeat (*). 

It turns out that the multiplication table T found by the above greedy 
algorithm yields another loop K of nilpotency class 3 whose inner mapping 
group is abelian. In addition, the following properties hold for K: (a) 
the nucleus is elementary abelian of order 16, (b) the left, middle, and right 
nuclei have order 64, (c) the two-element center coincides with the associator 
subloop. In particular, K is not isomorphic to L. Among other peculiar 
features, it contains a nonassociative power associative loop of order 64 that 
is the union of its nuclei. 

The construction of L takes a minute or so in GAP, since calculations 
in free groups are slow. A more direct, systematic, and much faster con¬ 
struction of L and K will be presented elsewhere [5]. 
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